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I. INTRODUCTION 
1.1 Background 

Although extensive ocean ambient noise studies were conducted 
during World War II with continuing effort since that date, few of 
either theoretical or experimental treatments have dealt with the 
directional properties of such noise. There are many possible sour- 
ces of ambient noise in the ocean any of which could produce non- 
uniform distributions. R. A. Prosch(!) in examining likely contrib- 
utors to the total acoustic energy content of the oceans has con- 
sidered, for example, biological sources, seismic noise, shipping, 
lightning and sea surface waves. His conclusion that sea surface 
waves are a major contributor is in agreement with World War II 
data. For example, the "Knudsen" (2) curves, a compilation of omni- 
directional measurements relating ambient noise spectrum levels to 
frequency as a function of sea state, exhibit a direct sea state/ 
ambient noise interdependence and by their very nature permit the 
inference of noise anisotropy in the ocean. 

Theoretical treatments considering the ocean surface as a 
possible noise source have been published by Eckart’), Roberson ^), 
Uriek(3) and Miller(9), Eckart has presented the phyeical picture of 
the ocean surface resolved into sinusoidal travelling waves of the 
particular phase velocity which would be associated with free flex- 
ural waves of an elastic plate excited by random local shocks of 
short duration. The above model would result in the generation of 
sound only dier the unlikely condition that the surface phase 


velocity exceeded the velocity of the bounded medium. Roberson 





without attention to a physical model of noise production considered 
both an infinite volume distribution of simple sources and the com- 
bined effect of a volume distribution of simple sources bounded by a 
surface distribution of dipole sources. Assuming an attenuation law 
proportional to frequency squared, he showed a -6 decibel/octave de- 
pendence of noise upon frequency for the unbounded case and a non- 
uniform dependence upon frequency for the mixed ease depending upon 
the ratio of the contributions from the two sources and the fre- 
quency. In this latter case, no combination of surface and volume 
sources produced a constant slope in agreement with the "Knudsen" 
curves except over a limited frequency range. 

Urick specifically considered the subject of directional 
distribution of ambient noise and concluded that significant aniso- 
tropy might exist depending upon the unimown manner in which each 
element of surface radiates. 

Miller considered the directional distribution resulting 
from the combination of straight path propagation from surface rad- 
lators close to the receiver and convergent zone propagation from 
distant surface radiators. The acoustic radiation from the surface 
varied as cosine Z, where Ô is the angle between the direction of 
radiation and a normal to the radiating surface. No bottom reflec- 
tion or volume sources were assumed. He concluded, tentatively, 
that noise in the above case is essentially isotropic over a hemis- 
phere with a slight increase within + 9 degrees of the horizontal 


and that the area immediately above and close to the receiver 





contributes about the same total noise as the contribution from the 
first convergence zone. 

If the ocean surface is a major contributor to ambient noise, 
it hag been reasoned that ambient level should decrease with depth. 
The change in ambient noise levels with depth 7) has been investi- 
gated experimentally. No such change has been detected to depths of 
700 feet in tests at 25 kilocycles per second. Urick!2) has shown, 
nowever, that the change to be expected between 3 and 300 feet depth, 
even at a higher frequency of 50 kilocycles per second where attenu- 
ation is greater, is only 3 decibels theoretically and thus a smaller 
change would have been obscured by measurement fluctuations at the 
lower frequencies, Teer (8) reported the results of tests at 19 
kilocycles per second with a small line array oriented horizontally 
and vertically in deep water. The results, which seemed to indicate 
that ambient noise is isotropic, are not considered conclusive be- 
cause of limitations in equipment and the environmental conditions 
experienced during the tests. 

In summary, no experimental data exist to describe the di- 
rectional distribution of ocean ambient noise. Theoretical treat- 
ments have indicated either isotropie or anisotropic distributions 
may exist depending upon initial assumptions as to the nature of the 


surface radiator and the uncertain effect of ocean sound channels. 


1.2 Research Objectives 


It is the purpose of this study to investigate anisotropy 


of deep ocean ambient noise in the frequency range of 750 to 1500 





cycles per second. Further specific objectives will be: 
a. To isolate vertical and azimuthal dependencies. 
b. To establish the identity of the sources contributing 
to non-uniform distributions, 


e. То infer the radiation patterns of such sources. 





II. AMBIENT NOISE EQUIPMENT 
2.1 Discussion 

There are two major reasons which have prevented more wide- 
spread investigation of the directional properties of ambient noise. 
The first concerns the limitation of a single directional hydrophone 
which, to monitor several directions in space successively, requires 
the assumption of a stationary proeess. The above difficulty may be 
overcome by the employment of more than one directional hydrophone 
or by forming several directional beams from an array of fixed non- 
directional hydrophone elements. Either solution generates a second 
problem, namely the difficulty of ensuring uniform gain in the sev- 
eral beams, Both limitations have been surmounted by a system oon- 
ceived by Anderson(?) and implemented by the staff of the Marine 
Physical Laboratory of the Seripps Institution of Oseanography, 
University of California, La Jolla. A 32-element array has been 
constructed in the shape of a great stellated icosohedron as shown 
in Fig. 1. A summary of the arrey design parameters is contained in 
paragraph 2.2 following. The 32 hydrophone outputs are combined to 
form 32 beams with axes coinciding with lines passing through the 
array center and the individual hydrophones, All 32 hydrophones are 
used in forming each beam. 

The problem of uniform sensitivity in each beam is minimized 
by employing clipped processing. At the output of each hydrophone 
channel prior to beam formation, all amplitude information is dis- 


carded by infinite clipping and phase information only retained. 
The clipped "on-off" form of the channel outputs permits digital 
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shift register techniques to be employed in the forming of the 32 
beams with an attendant saving in space and equipment complexity 
over conventional analog beam forming techniques. The system is 
thus independent of individual charmel gain differences from each 
hydrophone through the beam forming process. The problem of ampli- 
tude differences appears only in compensating for detector gain dif- 
ferences in the process of rectifying each beam output. The clipping 
and beam forming processes utilized in this system have been re- 
ported by Ander san‘ 10) end Rudniok 11), The restoration of the ana- 
log form of e signal from clipped hydrophone outputs is considered 


briefly in paragraph 2.5. 


2.2 Array Configuration Devel opment 
In arriving at the present array geometry of a great stel- 


lated icosohedron, atarsa) considered three factors: 

ae The number of elements, 

b. Element configuration. 

с. Element spacing. 
Since presumably small deviations from an isotropic field were to be 
measured, a symmetrical distribution was suggested permitting nearly 
identical directivity patterns of the several beams to be formed by 
electrical steering. A distribution of elements in depth rather 
than a spherical shell would permit greater interelement spacings 
for a given array diameter and thus increased array efficiency 
through minimizing interelement correlations. The number of elements 


wae limited by considerations of electronic complexity. Thirty-two 





appeared a practical number end (a) end (b) above were satisfied by 
the great stellated icosohedron. Optimum design of interelenent 
spacing was determined by calculations of array efficiency defined 
as the ratio of output to input signal to noise ratios. Signal in 
this case was the deviation from isotropy which it was desired to 
measure іп a given direction. An interelement spacing of 0.6 wave 
lengths at the 1 kilocycle per second nominal center frequency was 


chosen, 


Hw 
- 


2.3 Genera1 System Description 


An overall system block diagram is shown in Fig. 2. For 
simplicity only 1 of 32 identical receiving channels is represented 
completely. Reading from the left, an individual hydrophone output 
is first amplified and then employed to frequency modulate a tele- 
metering transmitter. The 32 transmitter carriers are distributed 
from 10 through 16.2 megaeycles per second in 200 kilocycles per 
second increments. An additional 1 megacyele per second amplitude 
modulated transmitter relays hydrophone array azimuthal orientation 
in accordance with the directional indication of a magnetic compass. 
The 32 hydrophones, preamplifiers, telemetering transmitters and 
campass system comprise the underwater portion of tho system. 

The transmitter outputs are coupled to the shipborne re- 
ceiving portion by eoaxial cable. A 1500 foot length was available. 
The 3/8 inch diameter, doubly armored eable serves also as the 
supporting member of the 700 pound receiving array. Cable attenua- 


tions of 1.5 and 1.8 decibels per 100 feet at 10 and 16 megacycles 
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FIGURE 2. AMBIENT NOISE EQUIPMENT BLOCK DIAGRAM 
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per second respectively necessitated the insertion of a broadband 10 
to 16,2 mececvele per second preemplifier prior to demodulation of 
the signal by the 32 matching frequency modulated receivers, At this 
stage it ie possible to view the amplitude and spacing of the 32 in- 
dividual earrfers as well as the frequency deviation imposed by each 
signal by means of a panoramic receiver and cathode ray tube indica- 
tor, 

Following demodulation in the FM receivers, the audio signal 
is further amplified and then infinitely clipped. А11 amplitude in- 
formation is discarded and only polarity information retained. The 
clipped gtznal then enters a digitel shift register ten bit delay 
line which provides suitable signal delay increments such that sub- 
sequent combination of the 32 signals in a mixing metrix yields 32 
beams symmetrically spaced over 47/2. The 32 beam outputs are then 
individually detected and integrated by a simple low pass RO filter 
of one second time constant. Following detection and integration, 
the 52 beam outputs are sequentially sampled at & 50 cycle per second 
rate by a mercury jet scanning switch for presentation on a cathode 
ray oscilloscope. Sixteen millimeter scope photographs preserve the 
data for analysis. 

All receiving equipment is housed in a portable Army Signal 
Corps communications hut. The entire 7,000 pound installation may 
be moved quickly by ecrane from its laboratory position to the fantail 
of PAOLINA T, a Seripps Institution of Oceanography research vessel. 
The electronic hut, array, associated winches and boom end gesoline 


powered 15 KW generator satisfying the 115 volt, 50 ampere demand 
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are shown in Figs. 3 and 4. 


2.4 Hydrophone and Transmitter Assembly 


The hydrophones are encapsulated barium titanate cylinders 
with a flat response averaging 90 decibels below 1 volt/y/ bar 
throughout the range of interest, 750-1500 cycles per second. The 
cylindrical sensitive element is 1-1/2 inches high, 1-7/16 inches in 
diameter and has a wall thicknese of 1/8 inch. Individual hydro- 
phone directivity is discussed in paragraph 3.5. Although the par- 
ticular units employed were constructed at the Marine Physical Lab- 
oratory, an identical unit, the REMACO Model R-130 is camercially 
available. They are rated at a working pressure in excess of 1600 
pounds per square inch, The hydrophone fits into the end of a 
16-1/2 inch long, 5/32 ineh thick stainless steel pipe which forme 
the watertight housing for preamplifier and transmitter. Ali hydro- 
phone units are series connected, 5O VDC required by each unit or 
1600 VDC total. A relay feature is included to remove a particular 
hydrophone unit should & vacuum tube filament open. Were this fea- 
ture not included the entire array would become disabled by the 
failure of a single tube, On the other hand, a short circuit any- 
where in the series etring disables all units beyond the point of 


shorting. A partially disassembled wit is shown in Fig. 5. 


2 Clipped Processin 
Because digital multibeam steering (DIMUS) is a relatively 
new technique, the following simplified illustration of the restora- 


tion of the analog form of a signal by the addition of several 








FIGURE 3. AMBIENT NOISE EQUIPMENT 
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FIGURE 5. HYDROPHONE / TRANSMITTER ASSEMBLY 
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samples of the signel which have been clipped is included. Referring 
to Fig. €a, the time varying analog audio output of an individual 
hydrophone channel is shown. Some band limiting has occurred in the 
hydrophone preamplifier bandpass filters. Further high-pass filter- 
ing by the simple network of Pig. 6b removes any remaining fluctua- 
tion below the frequency of interest resulting in the waveform of 
Fig. 6c. The waveform of 6c is transformed into one of two states 

as shown in Fig. 6d by the elipper amplifiers. Figure 7a repeats 
Fig. 6d. In addition Figs. 7b,-c, d and e represent the outputs of 
four other channels. The example has been so chosen as to illustrate 
the array gain of a five hydrophone system in detecting a directional 
signal mixed with omnidirectional noise. An input signal-to-noise 
ratio at the individual hydrophone of: unity has been assumed. The 
direct addition of the five output voltages produces a randomly 
fluctuating voltage as shown in Fig. 7f. The absence of a pronounced 
peak is to be expected since the noise at each hydrophone due to the 
isotropic background is uncorrelated while corretation between volt- 
ages generated by the directional noise will be visible only when 
the proper electrical delays corresponding to the geometrical delays 
between hydrophones are introduced by the shift registers. The 
voltages of Figs. 7a through с аге sampled at the shift register in- 
puts at a 5 kilocycle per second rate. This timing satisfies the 
Nyquist criterion of a minimum sampling rate equal to twice the 
highest frequency of the spectrum of interest. Figures 7g through 
7k repeat Figs. 78 through 7e with the proper delay introduced by 


the shift registers to correct for hydrophone geometry. The 
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FIGURE 7 THE ADDITION OF CLIPPED SIGNALS 





addition of Figs. /g through 7k shown in Fig. /1 again produces a 
randomly fluctuating voltage with the exception that now a peak ap- 
pears due to the directional noise. The reproduction of the original 
signal improves with the number of channels employed. In aural test 
with the 32 element system it is not possible to distinguish between 
an unclipped audio signal input as monitored at the output of a sin- 
gle hydrophone and the reconstituted input as heard at the output of 
a beam formed from the clipped signals of the 32 elements. 


^ 
- 


2.6 Array Beam Patterns 
The summed voltage output of all hydrophones may be written: 





|cos (%- G )cos(g-¢ }- cos(&- Z+%@ )cos(p- p+ е)! 


мһеге w += angular frequency, radians per second 


ғ,0,6 


12553 


hydrephone coordinates 


bes Ps = electrically steered direction 


a Po 
Due to hydrophone symmetry, a particular beam will exhibit conical 


angular deviations from steered direction 


symmetry about the beam axis. The icosohedron has two axes of вуш- 
metry and, accordingly, two different beam patterns. In a spherical 
coordinate system with the array center as origin, every element can 
be pictured ag positioned on one of two concentric shells. There are 
20 elements located on the outer shell and 12 elements on the inner 
shell. The two axes of symmetry are formed by lines through the 
array center and elements on the outer and inner shells. Since each 


array element has a diametrically located twin, the expression for 
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power output міћһ 6, = Ф, = О апі фо =0 Ъесалез: 


6 (9) [ez cos COS, 695 4 -cos(£- 6| Y Blarau 


where Q(w) ia a weighting function, the relative frequency response 
of the system and 4, 4, are the system bandpass limits. Computed 
patterns normalized to umit on-axis response are shown in Fig. 8. 
The continuous curve represents the solution to the preceding equa- 
tion for an outer beam (a beam whose axis passes through an outer 
element in the array). > 

The use of the preceding equation as written tacitly asaunes 
that the exact time delay specified by the expression =й созу COS Z 
is available. In practice the exact delays are not available owing 
to the limitations of the digital shift register which at its tapped 
outputs provides delays equal to integer multiples of the clock gen- 
erator period. The dashed curve in Fig. 8 is the outer beam re- 
sponse computed on the basis of the actual shift register delays. 
The dotted curve is the corresponding response for an inner bean. 
The combined effect of beam differences due to element configura- 
tions and steering delays results in an inner to outer on-axis beam 
response ratio of 0.853, a faetor which must be considered in the 


data reduetion process. 


2.7 Reduction of External Noise Interference 


The preliminary phases of at sea measurements confirmed the 
existanee of two sources of interfering noise which had been antici- 


pated. They were: 
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a. Shiv machinery noise transmitted through the hull struc- 

ture into the ocean. 

b, Self-noise generated by relative movement of the array 

and ocean water, 

Another troublesome source, external shipping interference, eould 
only be controlled by judicious choice of operating area free from 
shipping lanes or by testing during night time hours when shipping 
traffic wes at a minimum. A fourth source generated by the rolling 
motion of the testing ship is discussed in paragraph 5.1. 

Own ship machinery noi se was eliminated by two means. First, 
all normal ship mechinery was secured during the testing period. 
Secondly, interference from the special 15 KW gesolíne powered gen- 
erator serving as the power supply for the ambient noise equipment 
was reduced to a negligible value by isolating the generator from 
the ship deck by automobile tire casinge. The effectiveness of the 
isolation method was demonstrated by comparative evaluation with a 
special 5 KW sound isolated unit developed at the U, S. Navy Elec- 
tronics Laboratory for similar noise teste. This latter source was 
completely suspended from each of its four corners in a pipe support 
frame by airplane shock cord. Although its electrical capacity was 
not adequate for full time employment, alternate operation with 
first the limited capacity 5 KW wit and then the larger 15 KW 
source revealed no differences in ambient noise levels as measured 
at the array suspended in the ocean. 

Elimination of the second interfering source, self-noise 


generated by relative array/water motion, was accomplished by two 





means: (1) the suspension of the array from a system of buoye and 
floats in order to isolate the array from ocean surface and ship mo- 
tion and (2) the employment of a parachute sea anchor to eliminate 
ship drift caused by surface winds and attendant array/water rela- 
tive motion. The two techniques are shown schematically in Fig. 9. 
The suspension was во designed that the spar buoy would float in the 
position shown independent of array depth. The floats compensated 
for cable weight as the array was positioned at 100, 500 or 1000 foot 
depths. The spar buoy acts effectively as a spring with constant, 
K , equal to the buoyancy tame" ym length submerged. For the 4 inch 
diameter buoy in question / = 5 lbs/sec? which for an array dis- 
placement of 500 pounds results in a natural period of oscillation of 
[ = 2፦. /# 60 seconds 
a value 10 times greater than the longest period of the ocean waves 
encountered during the experiment. Particular attention was devoted 
to buoy swivels and connecting links in order to prevent noise gen- 
eration by any metal-to-metal contact. 

The sea anchor employed was a surplus government parachute 
fastened in the manner shown by manile hawser and equipped with a 
manila tripping line for collapsing purposes. Even in very light 
wind of 7 knots it wae effective in maintaining ship's head into the 
seas, reducing ship roll and completely eliminating drift. Under 
certain not infrequent conditions stretified oceen currents were ex- 
perienced where surface currents and wind were in a direction diffe- 


ent from the ocean current at, the array depth. The net effect wes the 
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FIGURE 9 SEA TEST CONFIGURATION 





towing of the array by the test shiv and parachute and invariably 


resulted in test cancellation. 








III. SYSTEM MEASUREMENT ACCURACY 
221 Discussion 

The ideal manner in which to specify system measurement ac- 
curacy would be in terms of beam response to a calibrated source cap- 
able of producing an output svectrum covering the ambient noise 
equipment bandwidth. Up to the present the above procedure has not 
been possible for a variety of reasons. The chief limitation arises 
from the low frequeney of operation of the system. Neither the nec- 
essary low frequency broadband source nor a quiet, deep water area 
with adequate facilities for ih array/source geometry has 
been available for test. Although it has not been possible to cali- 
brate the system as 2 single unit by one series cf measurements, 
calibration of all component parts of the system, individually, has 
been possible and is considered an adequate substitution for the 
single measurement procedure. In paragraphs 2.2 through 3.7 follow- 
ing, the possible error contributions of all component parts of the 
eystem are examined. Detector calibration has been established as 
the major source of error in the system and will accordingly be om- 
sidered first. It is the objective of the remaining paragraphs to 
show that all contributions from other components of the system are 


negligible in contrast to deteotor calibration orror. 


5.0 Detector Calibration 


The use of clipped processing considerably complicated de~ 
tector calibration procedures, Since the output amplitude of a 


given beam is a function of the relative phase correlation between 
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individual hydrophones and is independent of individual hydrophone 
output amplitude, the input/output characteristic of a given beam 
forming/detector combination could not be established by injeeting a 
signal of controlled amplitude into each clipper amplifier and mea- 
suring detector output voltage. 

The problem is further complicated by the nature of the shift 
register matrix output. <A representative waveform in shown in Fig. 
10a illustrating the summed voltage output from 42 shift registers 
with random noise input. The desired information is contained in 
the fluctuating AC. The — 5 kiloeyele per second clock 
pulses, arising from the particular form of shift register circuitry 
employed, eliminated the possibility of employing full wave rectifi- 
cation and necessitated use of less efficient half-wave rectifiers. 
A quantitative analysis of the system output including the effects 
of clipping is contained in Appendix I. 

Although 32 random noise generators whose relative correla- 
tions could be varied at will would permit detector calibration, no 
circuitry to fulfill this requirement was obvious. Instead a per- 
fectly correlated AO signal was employed. A 2.5 kilocyele per seomd 
test voltage of two phases, 180 degrees apart, was derived by fre- 
queney divider from and synchronized with the 5 kilooycle per second 
clock generator. Application of one or the other of the two phases 
to a particular clipper amplifier input would assure that the corre- 
sponding shift register outputs were identical, either off or on. 

By adjusting the number of clipper amplifier inputs of each phase, 


the matrix output for any beam could be varied from 0 to + 16 units. 
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Detector bias was adjusted so that with 17 inputs of one phase and 
15 of the opposite phase, the detector output amplitude was 1/2 the 
difference in clock pulse amplitudes as viewed by oscilloscope andas 
shown in Fig. 10b. Transient spikes as indicated in the figure and 
detector inefficiency prevented the achieving of exaetly zero output 
voltage for the condition of inputs divided equally as to phase. 

Time variations in individual detector characteristics were 
compensated for by trimming each detector bias to permit alignment 
with a predetermined test pattern as shown in Fig. i0c. This was 
created by steering the test signal to produce maximum response on 
beans 1 and 32. 

By progressively changing the inputs from half of each phase 
(minimum detector output) to all of one phase (maximum detector out- 
put) complete detector characteristics could be traced. Many trac- 
ings for each detector were required to average out the slight dif- 
ferences in individual register outputs emphasized by the artificial 
nature of the single frequency input but averaged automatically in 
practice by the random nature of the noise, The applicable portion 
of a typical curve is shown in Fig. 11. After one set of character- 
igties had been determined as described, а duplicate set was pro- 
duced. In по ease did the second set differ by more than 1 millivokt 
from the first set. This 1 millivolt error ig considered the maxi- 
mum system error and amounts to an error in the measurement of the 


output of any beam of 0.2 decibels. 
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4,5 The Атта 

Individual array elements comprising the 582302 receiving 
unit and the FM tranemitter canister housing were calibrated by the 
Us Se Navy Electronics Laboratory calibration facility at Sweetwater 
Lake, San Diego, California. Measured horizontal and vertical re- 
ceiving patterns were completely circular. Since the wavelength at 
the nominal array center frequency of 1 kilocycle per second is 5 
feet, the patterns are understandably uniform indicating that the 
individual hydrophones are equivalent to a point hydrophone. The 
cylindrical shape of the Huo Беше and canister does not alter the 
isotropic response. 

Although the response of the assembled array of individual 
elements could not be measured, acoustic theory and World War II 
measurement s (12) of sonar dome structures demonstrated that support 
members of cross-section equal to or less than 1/15 of a wavelength 
have no effect upon the beam patterns of arraye transmitting through 
such structures. Since the ambient noise array support members are 
1 inch steel angle iron approximately 1/8 inch thick, the above de=- 
sign eriterien is more tnan adequately satisfied. It is concluded, 
therefore, that the array may be considered formed of 32 exactly 
identical omnidirectional receiving elements unaffected acoustically 
by any supporting structure. Accordingly, it is considered that no 
errors are introduced by the array. 

3.4 The RF System 
The RF system was examined as to circuit noise, dynamic 


range and the introduction of undesired phase snifts. 
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With an appropriate capacitor substituted for the hydrophone 
crystal, electronice circuit noise levels averaged 20 decibels below 
sea state zero noise levels. 

The question of dynamic range primarily concerned chamel 
cross-modulation. Modulation levels were adjustéd such that maximun 
noise levels as measured at sea produced no more than a + 50 kilo- 
cycle per second deviation of the appropriate channel carrier. This 
adjustment permitted a 100 kilooycle per second separation between 
channels and eliminated all possibility of cross-talk between chan- 
nels. | 

Signal phase shifts wintentionally introduced by receiving 
channel components are equivalent to electrical time delays. Only 
differences between channel phase shifts are significant since equal 
phase shifts in each channel are equivalent to the application of a 
fixed time delay to all channels which would not affect beam output. 
Channel phase shift differences will not seriously affect beam out- 
put until they approximately equal the sample period of C.2 milli- 
seconds or 72 degrees. A large relative phase shift would cause the 
signal to appear at the incorrect shift register tap and would dis- 
tort the beam pattern. To measure the actual shift introduced by 
the system from hydrophone to shift register input, conventional 
Lissajous figure oscilloscope techniques were employed. An audio 
oscillator output was coupled by suitable attenuators to the oscillo- 
scope horizontal deflection plates. A parallel oscillator output 
was capacitatively coupled to the hydrophones in air by a special 


cup designed to fit concentrically over the array hydrophone 
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elements. The clipper amplifier output of the particular hydrophone 
channel undergoing measurement was fed to the vertical deflection 
plates of the oscilloscope. The oscillator frequency was adjusted 
for a straight line or zero phase shift. For the 32 channels the 
frequency of zero phase shift was found to vary between + 30 cycles 
per second centered at 1.01 kilocycles per second. A second meagure- 
ment of the phase shift characteristic of an individual channel 
determined by observing the frequency corresponding to the easily 
read Lissajous points of 180, 90 and zero degrees permitted conver- 
sion of the initial zero phase shift frequency measurement to phase 
shift at zero frequency. Fran this latter auxiliary measurement it 
was established that 0.33 degrees shift occurred per cycle or that 
the greatest measured shift was + 10 degrees from 1.01 kilocycles 
per second, a value well within the sampling period. 


It is concluded, therefore, that no significant errors were 


introduced by the RF systen. 


5.5 Beam Forming System 


Errors in beam output introduced by the beam forming system 
may occur due to four possible causes: 

a. Inadequate clipper amplifier dynamic range, 

b. Nonuniformity of clipper amplifier input band-limiting 

filters. 

с. Amplitude variation of shift register outputs. 


d. Nonuniformity of matrix summing resistors. 
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Clipper amplifier dynamic range, defined as that range where 
clipped output is independent of input amplitude, was measured to be 
70 decibels, a value well in excess of the range of ambient noise 
levels in the ocean. No error is attributed to this source. Ander- 
son(15) has shown that errors in output of individual elements com- 
prising an array degrade the array signal-to-noise output ratio in 
accordance with the relation: 

S/N = n(I1- o?) 
where 7 is the number of array elements and c is the standard devi- 
tion of the individual element ен The standard deviation for 
items (b), (c) and (d) above when measured individually and then 
combined was less than 10 per cent. Substitution of this value into 


the above expression ylelds a 0.05 decibels lose, a negliglble con- 


tribution to overall system error from beam forming components. 


5.6 Film Recording 


As noted in paragraph 2.3, individual beam outputs were 
photographically recorded continuously throughout each run by 16mm 
camera. Frame spacing wae geared down to 3.5 frames per second. The 
developed film was viewed for reading on a Craig 16m Projecto Edi- 
tor. The displayed frame size was as shom in Fig. 12 where a 40 
millivolt calibration pedestal, the base line and 32 beam outputs 
are indicated. An accurate, fine line, overlay grid of 4 millivolts 
per division was drawn to large scale and then photographically re- 
duced to the size of Fig. 12. The reading procedure first involved 


adjusting the projected frame size such that the calibration pulse 
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height corresponded exactly to 40 millivolts on the viewing grid. 
The value of each beam output could then be read quite simply to the 
closest millivolt. Assuming thet individual component errors add 
randomly, the error introduced in the film reading process is con- 
sidered small in comparison to the detector calibration error of 


paragraph 5.2, 


5.7 System Linearity 
A consequence of digital array phasing by the DIMUS technique 


and symmetry in the array beams is that the sum of the output powers 
of the 52 array beams ie a constant. Thus, an ineréase in level ina 
given beam will result in a corresponding decrease in level in the 
remaining beam outputs. A very strong signal, or its equivalent, 
good correlation in a given directicn, can result in beam capture 
with reshltant suppression of the remaining beam outputs and attend- 
ant loss of information as to the relative field intensities viewed 
by several beams. It is necessary, therefore, to determine that 
maximum output level on a given beam which does not introduce dis- 
tortion in the outputs of the remaining beams. Toward that end it 

is assumed that a signal-to®noise ratio at the input of an individual 
hydrophone of 1 to 10 represents a sufficiently small signal such as 
to preclude the possibility of signal capture. The output power gen- 
erated by such a signal-to-noise ratio has been derived in Appendix 
II and the results compared to maximum outputs observed during mea- 
surements at sea. Since no measured beam outputs exceeded those 


computed, it is concluded that the system behaved in a linear manner 








IV. TEST PROCEDURE 
4.1 Narrative 

The on-station test procedures were straightforward barring 
the occasional equipment component failure. Two hours prior to 
commencement of test the equipment was energized in order to permit 
a stable temperature to be reached. The detector circuits, for ex- 
ample, proved to be extremely sensitive to temperature changes 
greater than 1 degree and required thermocouple control of the two 
electronic equipment hut cooling fans. Upon reaching the test area 
the parachute sea anchor was rinsed and the ship permitted to reach 
an equilibrium condition headed into the seas. The array was then 
swung clear of the ship and lowered to within approximately 80 feet 
of the desired operating depth. The spar buoy and float system was 
then clamped to the array cable and the system lowered until the 
weight of the array was supported by the floats. The array was po- 
sitioned generally between 100 and 150 feet from the ship. 

Radar bearings and ranges to any contacts, enemometer read- 
ings and fathometer soundings were recorded prior to turning off all 
ship power with the exception of the ambient noise equipment genera- 
tor. The detectors were balanced to the calibration signal as des- 
cribed in paragraph 5.2 and an approximate 4 minute film recording 
made of background noise. Array compass heading, ship's head, swell 
direction, bearing to the array, oscilloscope and calibration pedes- 
tal and base line settings were recorded. 

Normally 3 recordings over an hour interval were made at 


each depth. At the conclusion of a given sequence of recordings at 
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one depth, ship's power equipment was energized, the array partially 
raised to permit detaching the float system and then lowered to the 
new depth. The procedure as described above was then repeated. Ap- 
proximately 6 hours were required to complete measurements at one 


station of 3 array depths. 


4.2 Test Environment 

Tests at sea were planned with the objective of maintaining 
a constant environment. Sinee ocean surface effects as noted in 
Section I were a likely source of non-uniform ambient noise distri- 
butions, an effort was made to control other possible variables such 
as the effects of land masses, surf, shallow water and varying types 
of ocean bottom. Although it would have been desirable with the 
above objective in view to eonduet all tests in the truly deep ocean 
beyond the continental shelf, only one such operation was undertaken 
owing to the limitations in ship speed which dictated a 24 hour one- 
way travel time to such areas. Only one such extended cruise was 
undertaken and accounted for all sea state three, 2000 fathom ocean 
depth data to be reported. The sea state one plus data were obtain- 
ed in 1000 fathom water in an area which shoaled to 500 fathoms in 
all directions within approximately 10 miles. The sea state one- 
half data were obtained in a 1000 fathom area similar to that 
studied under sea state one plus. 

Sea state classifications are at best an approximation even 
when derived from a stable platform. Since the difficulty of main- 


taining an objective viewpoint from a rolling and pitehing 75 foot 
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test platform was appreciated, the initially recorded sea state 
values of two and four were reduced to one plus and three respec- 
tively following examination of motion pictures taken at the time of 
recording. All observations were conducted in moderate swell con- 


ditions. The environmental conditions are summarized as follows: 


Sea State Ocean Depth Array Depth Type Bottom 


2 2000 fathoms 150, 560, Gray mud 
1000 feet 
1+ 1000 fathoms + 150 feet Grey mud 


1/2 1000 fathoms 1000 feet Gray mud 





V. MEASURED DATA ANALYSIS AND RESULTS 
5.1 Noise Source Origin 

Non-isotropic noise distributions, functions of both azimuth 
and elevation angle, were measured during the tests deseribed in 
paragraph 4,2. Observations of the charactor of the noise fluctua- 
tions in the upper beams of the array, namely the vertical upward 
beam directed normal to the ocean surface and the next family of 5 
beams at 35 degrees from the vertical, indicate that one of the sev- 
eral noise generating mechani ems is the white caps of bresking 
waves. Near the ocean surface individual beam outputs vary consid- 
erably from one another in a random manner, Momentary increases in 
level in a given beam were probably due to the arrival of energy 
from en individual or small group of white caps or wavelets. As 
array depth increased, beam-to-beam differences decreased. The 
logical explanation lies in the increased area of ocean surface sub- 
tended by the beams and the liklihood that the beams of a particular 
elevation angle view more nearly the same distribution of white caps, 
The above conelusion was substantiated during tests designed to re- 
dues external noise interference described in paragraph 2.6. 

Despite efforts to sound isolate the 15 KW power generating 
equipment during early phases of the experiment, interference from 
the test ship invariably was observed in that beam most nearly di- 
rected at the ship. When installation of special sound isolated 
generating equipment failed to reduce the interferenee, further ex- 


amination revealed that the teet ship was unwittingly generating the 
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equivalent of white caps. The heavy handling boom, rigged to port 
for all tests created a pronounced list. The list, coupled with a 1 
to 2 foot freeboard amidships, resulted in the deck edge dipping be- 
low the waves as the ship rolled even in the lowest sea states. A 
considerable volume of water entered through deck ейде seuppers. On 
the return roll this water would pour from the seuppers on to the 
ocean surfece like water from a pitcher. The process was repeated 
with each roll. The interference wes substantially reduced and the 
source identified thereby by streemiining the entire portside, 
rounding over projections and soveriné scuppers with a 50 foot long, 
4 foot wide strip of canvas. 

White caps are probably the principal eontributor to back- 
ground noise. However, the ocean surface by some unexplsined mech- 
anism still contributes to the noise field at sea states above zero 


even when white caps ere not present. This fact is substantiated by 


the measured distributions reported in paragraph 6.3. 


T Data Analysis 


The film data was first reduced to tabular form as described 
in paragraph 3.6. Three equally spaced 4 minute samples recorded 
within a 1 hour interval and with an observation every 14 seconds 
within each sample resulted. The 3 samples for each array position 
were examined for fluctuations within each sample and between sam- 
ples. Fluctuations within samples varied from less than + 0.4 deci-~ 
bels for beams directed in the lower hemiephere to + 1.0 decibels 


for beams directed upward. The standard deviation of sample 
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averages was smaller generally than the standard deviation of obser- 
vations within a sample by the factor = ша where 7 is the 
number of observations oomprising a sample. Accordingly three sem- 


ples at each depth were averaged and treated as one. 


5.5 Azimuthal Dependence 


The data were first analyzed for azimuthal dependence by ex- 
emining the outputs of families of beams of a particular elevation 
angle. In general, some azimuthal dependence was apparent princi- 
pally in the outputs of those array beams directed in the upper 
hemisphere. The distributions normally were elliptical in shave. 

A representative azimuthal distribution is shown in Fig. 15. It is 
considered simmificant that the maximum beam outout for a given elə- 
vation angle normally aligned itself with the wave fronts of the 
swells while a outputs usually were pervendicular to ths 
swells, This fact suggests a non-uniform radiation pattern of the 
white cap or wavelet source; transmission іп directions normal to 
the wave front 1s impeded by the trouchs, with lateral transmission 


along the ridges being less affected, 


5.4 Vertical Dependence 


The data were examined for non-uniformity in the vertical 
plane by averaging the outputs of all beams of a given elevation 
angle, thus eliminating azimuthal dependence. The results are shown 
in Fig. 14 which presents the measured field as a function of eleva- 


tion angle for the conditions listed in paragraph 4.2. Sea state 
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FIGURE 13 AMBIENT NOISE AZIMUTHAL DEPENDENCE 
FOR VARIOUS ELEVATION ANGLES(8) 
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FIGURE 14 VERTICAL DISTRIBUTION OF AMBIENT NOISE 
VS SEA STATE AND ARRAY DEPTH 





effects are apparent with meximum to minimum differences for a given 
condition varying from less than 1 decibel for sea state one-half to 
б decibels for the highest sea state. Less simificance can be at- 
tributed to differences in anisotropy magnitude due the variable of 
array depth because of uncertainty as to the constancy of sea state 
during the recording period. Approximately 6 hours were required to 
record data at the three depths during en early morning period of 
veriable winds. It is not unlikely that the recorded sea state of 
three could in fact have fluctuated between sea state two plus and 


ቁብ. 
> 


three plus. 





VI. TRUE FIELD RESTORATION 
6.1 True Field Indeterminacy 
In the preceding paracravhs the analysis has been confined 
to the measured fields. It is necessary to distinguish between the 
measured distributions and the true noise distributions. The output 
from a beam whose axis is oriented in a fixed direction is not de» 
termined solely by the true field intensity in the oriented direc- 
tion. The beam output is propertional to the weighted mean of the 
intensities in all direetions. The weighting factor is the beam 
pattern, Mathematically expressed, the output of a beam becomes: 
Y [ 614-6, 6 $|x ай (1) 
where Y = beam output or measured field. 
G = the beam response, a function of the angular difference 
between axis orientation, 0 , Pg» and a direction of 
interest 2, 2 
X = the true field intensity where, for simplicity, azimuthal 
symmetry in the noise distribution is assumed. 
Since the array beam patterns are not pencil-like with zero response 
in all but the desired direction but are broad and possess a side- 
lobe structure as shown in Fig. 8, the effect is to smeer the fea- 
tures of the true field. 
The problem of reconstructing the true from the measured 
distribution is common to the field of astronomy where a radio an- 
tenna is used to survey the distribution of radio brightness over the 


sky. The arguments of pria who first considered the problem 


ኣሪ 
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in detail in conjunetion with aerial smoothing in radio astronomy may 
be adapted to the underwater acoustics case as follows. Upon inte- 


grating over the variable, 4 , equation (1) becomes 


7/2 

Y (4) Ди E N (2) 
where G signifies the average of G over all azimuth for a given 
colatitude and X'= X (@)sınd. The integral on the right, a con- 
volution integral, may be written 

y (s) = PREX (2) 

By means of the convolution theorem, equation (3) may be rewritten in 
terms of the Fourier transforms of the three variables as follows 
Y(s) = G(s) X(s) (4) 
where the convolution symbol has been replaced by the algebraic pro- 
duct and the three terms are the spectrums of the corresponding func- 
tions of equation (3) and "S" 4s the number of cycles per unit of 
angle, 2 . From (4) it is seen that У (5)=0 whenever С (5)= О 
and thus X (s) is indeterminate under such conditions. Obviously 
in any restoration process dependent upon measured values of Y only 
the low frequency components of X are restored to their full value. 
The resulting values of X have been called the "principal solution" 
which may be defined as that solution whose spectrum is the same as 
the spectrum of the true distribution at all values of S for which 
G(s) ЖО, The results of a restoration process to be described in 
the paragraphs following, are shown in Fig. 15. The significant 
feature of Fig. 15 is that the trensform of the restored field 
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approaches zero at almost the same rate as the transform of the mea- 
sured field, This is an important condition since any values of \(s) 
for S greater than S, , the value of 5 where G(s)= O, would rep- 
resent fictitious invisible distributions introduced by the restora- 
tion process. Since ZX (s;  О , the addition of fictitious high 
frequency eomponents is negligible and X (Q) may be considered an 


exact principle solution throughout this study. 


6.2 The Matrix Approximation 


The restoration process employed ean be, considered in two 
parts: the estimation of a first approximation, and an iterative 
solution deriving from such first approximation. The former, the 
matrix approximation, is novel and is considered in this paragraph. 
The latter, similar to methods of Bracewell, is described in para- 
graph 6.5. 

The matrix approximation is generated by converting the con- 
volution integral of paragraph 6.1 to a finite summation of terms of 
number equal to the number of beams. Each term in the summation rep- 
resents a contribution from a particular solid angle and is propor- 
tional to the average of the true field subtended by the solid angle. 
Since a similar expression can be written for each beam, a square 
matrix results relating values of measured field to the unimown 
values of the true field. By inverting the matrix, the true field 
values can be obtained. An exact matrix can be written only for the 
isotropic ease. However, the matrix so derived when used with поп- 


isotropic distributions yields true field values which although 





inexact, are closer epproximations to the true field distributions 
than ean be arrived at by intuition. The approximation so derived 
may be then conveniently operated upon by iterative techniques to 

yield the principal solution. The details of the derivation of the 


matrix appromimation are contained in Appendix III. 


6.3 The Iterative Solution 


The iterative solution is a method of restoring the true 
field by successive substitutions. The first approximation pre- 
viously described is tested by numerically completing the integra- 
tion of equation (1) of paragraph 6.1, and comparing the results with 
the measured values of beam power, Y, +. From experience with num- 
erous cauputations of this nature, it has been determined that a 
suitable second approximation is obtained by correcting the first 
approximation by an amount equal to twice the difference in decibels 
between the measured and computed values of Y . Three such itera- 
tions were usually sufficient to reduce the difference between mea- 
sured and computed values to less than experimental measurement 
error. 

Bracewel] (14) has shown that a sufficient condition cf con- 
vergence of the above process it that G(s) be non-negative, a condi- 
tion which is satisfied by the transforms of the bean patterns of the 
experimental array. 

A comparison of the decibel differences between a given pair 
of beams before and after restoration show a two to three-fold in- 


crease in the decibel difference. It is reasoned that measurement 





errors will be similarly increased in the restoration proeess. Ac- 
cordingly, the errors in the true field principal solution are ев- 
timated to be no greater than + 0.6 decibels except in those regions 
where abrupt discontinuities in the distributions exist. A sample 
computation of the complete restoration procesa is contained in 
Appendix IV. 

The resulting true field distributions for the five environ- 
mental conditions described in paragraph 4.2 are shown and tabulated 
in Figs. 16 through 20. The existance of a marked noise null in the 
horizontal direction is ሀ9786866:11.ቀ useful in the design of sonar 
arraya whose objective is to detect weak signals in a noise back- 
ground. For the detection of horizontally arriving signals in a 
background of noise arriving from all directions, good vertical 
directivity as well as horizontal directivity is advisable. To 
illustrate: consider two arrays of equal horizontal directivity but 
differing vertical directivity. Let one array have uniform vertical 
response and the second 20 decibel rejection over all vertical anges 
except within 15 degrees of the horizontal. The signal-to-noise 
ratio of the latter array will be 5.5 and 8 decibels better than the 
former when both sre positioned in an isotropie noise field and the 
field of Fig. 16 respectively. 

As previously noted, variations in the noise field distri- 
bution for differing array depths could not be established because 
of uncertainty as to the constaney of sea state during the recording 
period. The increase in anisotropy with increasing sea state, how- 


ever, is obvious. 
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FIGURE 18 TRUE FIELD VERTICAL DISTRIBUTION 
SEA STATE 3 ARRAY DEPTH 150 FT. 
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FIGURE 20 TRUE FIELD VERTICAL DISTRIBUTION 
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VII. AN OCEAN MODEL 
1。1 Ocean Model Development 

From the nature of the distributions shom in Figs. 16 
through 20 and consistent therewith, the following simple picture of 
the mechanism of noise generation and redistribution in the ocean 
has been constructed. Noise in the deep ocean is presumed to be a 
superposition of two fields. The first field ia that which exists at 
sea state zero prior to the development of wave motion and/or white 
caps. Its distribution is isotropic and its origin was not estab- 
lished during this experiment, а жагыз» explanation of the origin 
of low sea state ambient noise is considered in paragraph 8.1. The 
second source of noise is generated by wave motion and white caps at 
the ocean surface. The equivalent surface radiator is highly direc- 
tional in the vertical plane when azimuthal effects are averaged. 
The directionality can possibly be attributed to the ocean waves and 
swells functioning in a manner analogous to engular focusing reflee- 
tors. In addition, in azimuth the ocean swells impede transmission 
at shallow transmission angles in directions normal to the wave 
fronts of the swells. 

It is further postulated that the noise intensity at a given 
point in the ocean due to the equivalent surface radiator arrives by 
both a direct ray path and single bounce path experiencing specular 
reflection at the ocean bottom. The intensity at the surface source 
experiences inverse square spreading loss and absorptive attenuation 
of the tem 10 የ2 where 2 28 ቂከፀ absorption coefficient 


(decibels/kiloyard) and / is ray path length in kiloyards. 


27 
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Referring to the model of Fig. 21 and disregarding fer the moment 
the low sea state isotropic contribution to the total noise field, 
the above postulation of an ocean model may be restated in methe- 
matical form as follows: 


-0 fp 


X (8) dN = E 0 LN (1) 


D 


where Х (0) 012 = noise intensity per incremental solid angle 
at the receiver in the upper hemisphere 
S(@)dA = noise intensity at one yard per incremental 


® 
> 


radiating surface. 





Since 
40 - 222282 (2) 
D 
equation (1) may be rewritten es 
CAG 
КО - SX s то 19 (3) 
cose 
or, finally, 
q^ sec 
S(0)= X(9) cos 9 10 0 О<й< 7/2 (4) 


A similar expression for the intensity at the receiver in the lower 
hemisphere due to bottom reflection may be written 
AMF AN 


90) 2да (6) |0 О (5) 


pecus gf. 
aa Е аа 


where z (2) is the intensity bottom reflection coefficient and 


is the supplement of Fo e 
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FIGURE 2! OCEAN MODEL L 








Since 
еке d Atte, (6) 
(Mp + Aro) 
equation (5) may be written as 
ale D-H) 52” 
a(6@)= да cos 9 10 ы 6 ааа а. (7) 


Да < 7 


1.2 Dedueed S (9) and a (9) 


Employing the distributions of Figs. 16 through 19 and equa- 


tions (4) and (7) above, values of 5 (2) ап4 с (2) меге computed. It 


was necessary to estimate the contribution of the isotropie zero sea 
state noise field and to subtract that amount from the total field, 
X (0), prior to proceeding with the computation, Several alterna- 
tive methods of estimating the zero sea state contribution, none of 
whieh can be considered exact, are possible. For example, the level 
of the null oecurring at an elevation angle of approximately 100 de- 
grees in all distributions may be considered due to the low sea 
state isotropic field. Errors in this assumption are likely due to 
the difficulty in establishing the true minimum. The most probable 
point of error in the restored distribution is at the minimum due to 
the difficulty inherent in resolving a sharp discontinuity with e 
broad beam pattern array. In addition, contrary to the very simple 
straight ray assumptions of the ocean model, curved path arrivals 
direct from the curface or by convergent zone paths 5) due to verti- 


cal velocity gradients probably contribute to the null level. For 
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exemple, 8611(12) has considered the deep ocean case with negative 
gradient to deep depth. He reports bottom reflection losses for 
grazing angles of incidence which are approximetely 6 decibels less 
due to the refractive effect of the negative gradient than would be 
predicted on the basis of straight ray path assumptions. 

The best estimate of the isotropic contribution utilized the 
Knud sen ^) sea state/noise level curves. If the restored distribu- 
tion is averaged over 4r 2 , the relative level derived from an 
omnidirectional hydrophone can be obtained. Subtraeting from this 
level the decibel difference TUYN the Knudsen noise level for the 
sea state under consideration and sea state zero, the level of the 
isotronte eontribution can be established. The weakness in the above 
approach rests in its dependence on subjective estimates of sea state 
level. The error introduced by this technique decreases with 1п- 
creasing sea state since the percentage contribution to the total 
field of zero sea state isotropic noise becomes smaller under such 
conditions. 

Figure 22 18 a plot of S(f), the surface source radiating 
pattern, deduced by the above procedures from the restored fields of 
Figs. 16 through 20. A reasonable approximation to5 (/)is the 


fimetion: 
90-2 


аа 


5(0) - (cos 4) 19 O PO 
S(9)= cosé 9 70° < @ < 90° 


which 4s shown superimposed on the experimental data. 
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The bottom refleetion coefficient (reflected pressure/ 
incident pressure) derived for conditions corresponding to Fig. 22 
are plotted in Fig. 23 together with experimental data and theoreti- 
cal plots from MacKenzie! '6) and Marsh 17). The fit is eonsidered 
good in view of the indirect manner in which the data were obtained 
and the several assumptions involved in the analysis. The results 


are considered to substantiate the simple ocean model proposed. 
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VIII. SUMMARY 
8.1 Conclusions end Diseussion 

Prior to this experiment, as noted in Chapter I, no experi- 
mental data existed to describe the directional distribution of deep 
ocean ambient noise. Accordingly, stated objectives of this study 
were to investigate anisotropy of ocean anbient noise, to isolate 
vertical and azimuthal dependencies, to establish the identity of 
sources contributing to non-uniform distributions and to infer the 
radiation patterns of such sources. All the above research objec- 
tives have been realized by this experiment. Deep ocean ambient 
noise is anisotropic in the frequency bend 750 to 1500 cycles per 
second and at sea states above zero. Vertical and agimuthal depen- 
dencies have been isolated. The surface of the oeean has been 
identified as a major contributor to non-uniform distributions and 
the radiation patterns of such surface sources have been deduced. 

Vertical anisotropy arises from the superposition of a 
noise field due to the radiating ocean surface upon an isotropic 
field present in the absence of surface radiation. The shape of the 
vertical distribution, Figs, 16 through 20, may be compared to the 
form of an arrowhead. The arrow tip is directed normal to the ocean 
surface and represents maximum intensity arriving from the surface 
direction. The arrowhead tail indioates noise originating at the 
ocean surface which has been specularly reflected from the ocean 
floor. A dip in the arrowhead tail in a direction normal to the 
ocean bottom is evidence of a lower value of bottom reflection co- 


efficient for normally incident sound than for other angles. A 
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marked null at the waist of the arrowhead, approximately contained 
within + 15 degrees of the horizontal, exists in the vertieal dis- 
tribution. The mamitude of the null varies between 2 and 12 deci- 
bels below the maximm intensity normal to the surface for sea states 
one-half to three respectively, The null is due to two effects. 
First, contributions from the surface source to the horizontal or 
near horizontal distribution arrive from greater distances than 
contributions to the off-horizontal distributions and have exper- 
ienced, therefore, greater attenusetion fron spreading and absorp- 
tion. Second, the level of — fran the surface source in 
horizontal directione ia more than 20 decibels below the level of 
radiation in the downward directions. 

The non-uniform noise distribution with a minimus in the 
horizontal direction is of practical importance to the design of 
sonar receiving arrays. The importance of vertical directivity to 
the detection of horizontaily arriving signals is apparent. 

The origin of the isotropic zero sea state noise field was 
not established by this experiment. Of the various sources con- 
sidered by Frosch V!) distant surface motion is considered most 
likely. Although the sea surface may be flat at the point of mea- 
surement, the noise produced by distant surface disturbances can 
propagate very long distances with negligible attenuation due to 
absorption in the 750 to 1500 cycle per second frequency band ob- 
served. The scattering and reflection processes experienced by the 


gomd in reaehing the point of observation could destroy its origi- 


nal charactor and an isotropic field result. 





Ocean swells are probably responsible for azimuthal aniso- 
tropy br impeding transmission of surface generated noise in dirce- 
tions normal to the weave front of the swells. 

The intensity at a receiver due to a unit incremental sur- 
face radiator of intensity,S(@), at one yerd, can be considered ag 
the sim of intensities due to direct transmission and specular bot- 
vom reflection with attenuation due to spherical spreading and ab- 
sorption occurring over the separation between radiator and receiver, 
The quentity, S( 7), ic markedly directional in the downward di- 


® 


rection and can be represented approximately by the expression: 


20-0 
/ሪ 


(cos @) О< д < 70° 


S(O) = 
cos* 2 709 « д< 90° 


The above ocean model is not unique, It does, however, con- 
form to the measured data as evidenced by the reasonable agreement 
between bottom reflection coefficients determined from this experi- 
ment and values reported by other investigetors. Further additions 
to the model, as for example the consideration of bottom/surface 
multiple reflections, are obvious. Such refinements would presume, 
however, a knowledze of bottom topegraphy and array position rela- 
tive to the bottom greater than existed during the experiment. 

Although the data upon which the ocean model is based were 
obtained in ocean depthe greater than 1000 fathons, there is no 
known feature of the model which would preclude its application to 


shallower depths of perhaps 100 fathoms. Ite application to oven 





shallower depths will depend upon the presently unknown effeet of 
water depth upon the nature of the surface source radiation pattern. 
On the basis of measurements conducted over a one hour per- 
iod under constent conditions of array orientation and sea state, 
the ambient noise field can be considered stationary. A single, 
trainable and tiltable directional receiver could be employed, 
therefore, to duplicate the performance of the 32 element, nulti- 
beam array. However, to achieve beam patterns equivalent to those 
of the 32 element array with a single directional unit, focusing 
reflectors of dimensions А to the 32 element array would 
be required. А very difficult problem of mechanical training and 
elevating the reflecting unit would be substituted for the compara- 
tively simple telemetering problem associated with the 42 element 


array. 


8.2 Suggestions for Further Experiment 


The results contained herein represent only a first step in 
the area of ambient noise directional distribution investigation. 
Objectives for future experiment should include: 

ae Shallow water investigation. 

be Extension to higher sea states. 

с, Extension to array depths of at least 5000 feet. 

d. Extension of bandwidth covered from less than 400 

cycles per second to 5 kilocycles per second. 

e. Investigation of the effect of land masses such as 


islands and the continental shelf and the influence of surf. 
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An additional area of investigation is to substantiate the- 
oretically and experimentally the rediation pettern deduced from thia 


experiment of a white cap situated on the crest of a swell. 
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APrENDIX I 


QUANTITATIVE SYSTEM ANALYSIS 


The objective of this section is to derive an expression for 
the voltage appearing at the output of a given beam deteetor for a 
random noise input. 

Consider first an /V element array of omnidirectional elements 
situated in an isotropic noise field of purely continuous spectrum. 
If the hydrophone electrical outputs generated by the noise field are 
summed, with or without the introduction of an electrical delay, and 
if the hydrophones are sufficiently widely spaced such that cross- 
correlations between elements are zero, the array output following 
square law detection will be proportionel to /V к . If the eross- 
correlations are not equal to zero, the output will be proportional 
wa (1 Awe ye where 


| N(N-1) 
Өт Дои 
IF) Y 


and P, is the normalized cross-correlation function derived from 
the input noise spectrum. The value of е “may also be determined 
from the directivity ratio. Utilizing the beam patterns of paragraph 
2.6, the directivity ratio is derived by numerical integration from 
the expression 
pè ~~ ._ 
| 609, 6) 2.0 


where G (0, Ф) ів the beam pattern of paragraph 2.6, and G(0,0)= |. 


10 





ያ! 


The quantity, о , шеу be determined from the expression 
Wit @=|/OFR. For inner and outer beams c" has been com- 
puted to be 0.25 and 0.02 respectively. Lawson 17) has shown that 
the correlation functions of unclipped and clipped signals are re- 
lated Ex A, 


clippe 
as is the case under consideration, the arcsine of the argument is 


¿2/7 arcsin (ae ae For small values of correlation 
equal to the argument and the correlation function of the clipped 
signal reduces to 2/7 times the unclipped correlation function. 
Since the actual detector employed was a linear rectifier rather than 
a square law device, an added factor K » 0.8 is required. K=J 9 
is the ratio of the mean absolute value to the root-mean-square 
value in a normal distribution. A further factor of 0.5 is intro- 
duced because the scale of the abscissa of Fig. 11 is 16 for maximum 
detector output rather than 32. The detector output voltage is 
therefore proportional to 

(2) (xe) 


Substitution of the values V = 32 and е = 0.25 into the above ex- 


маг 


pressions produces the value 2.44, This quantity when entered as 
the ebscissa of Fig. 11 indicates an output value of 52 millivolts 
for detector No. 1, an inner bean. 

The computed value of 52 millivolts was experimentally con- 
firmed from data recorded at sea. When the outputs of all 32 beams 
are averaged, the resulting value is equivalent to an individual 
beam output to be expected from an isotropic, random noise field. 
The averages of several sets of data equalled approximately 52 milli 


volts. 





APPENDIX II 


SYSTEM LINEARITY 


Rudnt ek 11) includes the necessary expressions to compute 
beam output power ns a function of array geometry and input signal- 
to-nolse ratio at an individual hydrophone. Output power may be 


expressed ав 


Me 
n Me 
ጌ 


~ 
э 
1 
һ 
чм 


(1) 


with the summations to be completed over the /V array elements. The 


quantity, A (O), is the combined autocorrelation function of 


„ 
signal and noise. For purposes of this analysis, the spectra of 
signal-to-noise may be considered identical and are assigned dif- 
ferent nomenclature only to differentiate between that noise which is 
within the main lobe of the directed beam and that whieh is outside 
the lobe. Implicit in the last statement is the assumption that 
noise is completely uncorrelated with respect to angle of arrival. 

со „ (O) is the autocorrelation function of the spectrum introduced 
by clipping and is defined by Rudni ck ! !) as 


2 (m'a, c1] E фар (е) е, (г) (2) 


^ (O) may be expanded to 


/7 
O = С 
MAS г. P .5) d 








257 (О) 5 
(3) 


me 





15 


where o and o are noise and signal autocorrelation functions re- 
spectively, с, (О) = | since by definition it is the signal to 


ሃ ኘሣ 
which the beam is steered. У 7 c, (O may also be separated into 
/ =| yl š 


N N MN“) 

> o ER 

E cu - 2 QUEE 22 €, (0) (4) 
It follows fron (2) that С = 7/2 -|. In addition, Rudnick has 


shown that 
// 


ቿጅ ር O) 
'ፖ/ А 


is small and may be neglected. 
Substituting the several expressions described above into 
(1) and reducing terms, the output voltage of the linear half-wave 


rectifier can be shown to be proportional to 
1/2 да 
በቤ. -] ا‎ aoe (//+1)| (2) 


If a small signal-to-noise ratio input of S/n equal to 0.1 
is assumed and with M egual to 52 and "e" equal to 0.25, equation 
(5) yields the value 3.87. Entering the detector characteristic of 
Fig, 11 with the above value, an output voltage of 90 millivolts is 
obtained. No recorded data was used in the analysis where the out- 
put of any beam exceeded 90 millivolts. It is concluded, therefore, 
that the system was operated throughout its linear region during the 


tests. 





APPENDIX III 


DERIVATION OF THE MATRIX APPROXIMATION 


Considering the case of azimuthal symmetry, the output of 


any vertical beam may be written: 
¥ =f x8) G, (8,9) 20 (1) 
477 


where Y «= output of ,” beam 
X = continuous true field distribution (power/unit 
receiving area/unit solid angle subtended by the 
source ) 
G = power response of »” beam 
If the true field is assumed to be isotropic, then 


y = X, | G, (8, $148 (2) 


77 


The directivity ratio of an array is defined as: 


| Ј 6,02, é) 40 


DR (3) 
G, (2, A) | dN 
47 
Combining equations (2) and (3) 
y = Xo G ао ат (4) 


where G (A, Ф) equals unity and 0.853 for outer and inner beams 
respectively as defined by paragraph 2.6. Setting 47 Xp = | 
which has the physical interpretation of total power from all direc- 


tions normalized to unity, then 
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D.R. (outer beam) 
- (5) 
0,855 D.R. (inner beam) 
Rewriting equation (1) as a sum of eight terms, with each term rep- 
resenting the contribution from a spherical zone formed by two 


planes parallel to the equator and each zone containing one of the 


eight vertical beams, then: 


7-8 27 
| ak J sin #22 ] E (6) 


ገሙ 


~ 
= 


The true field, X, , may be non-isotropic, although ላ) is a constant 
within a particular zone. The quantity G, (0, $) may next be aver- 


aged over all azimuth. 


731! 


J-B 
g” 27 6 (9 2(0) sinddé (7) 
ё 


where @ апд 2 аге the limits of the spherical zones. There 


exists between the limits of each zone a particular value оР 0:70 


which permits equation (7) to be transformed to the following 


7 


ао С 28 24% (8) 
| 


/:8 
J=1 
Equation (8) may be further simplified by replacing the integral by 
a constant, К, ‚ for a particular zone yielding 
8 


/ 
y = 27 С ЦЕНИ 2 K X (9) 


VES 
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or more briefly 


Y, = Фу Х, ly А = (8 (10) 


Equation (10) relates the measured field values, Y, , toa set of 
eight simultaneous equations in eight unknowns which are the true 
field values, X, . 

The 0, matrix which satisfies equations (5) and (9) was 
determined by trial and error. The matrix was then inverted by 
digital computer to produce a set of equations of the form 

X= a, Y, 
The coefficients of the inverted matrix are listed below, By mul- 
tiplying each element of row one by Y, , row two by Y> , etc., the 
sum of all oolumn one elements ylelds X , the column two elements 
X>, tc. The values of X so derived are the first approximation 
to the true field in the restoration process. 
Table of Coefficients a,, 


// 1 2 5 А 5 6 T 8 
1 67.855 -32.537 27.269 -19.912 10.678 -19.322 5,122 -8,932 


2 -1.677 89.327 -67.679 23.059 -19.345 16.734 -5.158 -1.845 
5 20.056 -9.622 92.536 -39.901 20.905 -32,688 0.808 -0.282 
h «1,102 2,505 -52.047 106.681 «38.812 20.422 -5,628 -0.799 
5 0.799 -5.628 20,422 -38.812 106.681 -52.047 2.505 -1.102 
6 -0.282 0.808 -32.688 20.905 -39.901 92.536 -9.622 -0.056 
7 1. B85 «158 16.734 -19.345 23.059 -67.679 89,327 -1.677 
8 «8,932 5,122 -19.322 10.678 -19.912 27.269 «32.537 67.855 
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APPENDIX IV 


TRUE FIELD RESTORATION SAMPLE COMPUTATION 





In order to illustrate the restoration process, the proce- 
dures followed in deriving Fig. 18 will be described. 

Table I is a random selection of recorded data typieal of 
that which results following the film reading procedures of para- 
graph 3.6. The values shown are the individual beam detector out- 
puts in millivolts following one second integration end at 14 second 
intervals. Beam numbering is such thet beam 32 is directed at the 
surface, beams 27-%1 at 36 degrees from the vertical, beams 22-26 
at 64 degrees from the vertical, beams 17-21 at 78,5 dezrees from 
the vertical, etc., to beam 1 directed at the bottom of the ocean. 
The RMS values for each beam for each rum are computed and tabulated 
in columa 2 (Aver) of Table II, The detector correction is applied 
(Corr Aver) and the RMS output for beams of a given colatitude сош- 
puted (Vert Aver). The average of all beams ie similarly determined. 

Power ratios fer each colatitude with respect to beam 52 and 
with respect to the average of all beams are computed. The Y 
column is determined by multiplying the appropriate expression of 
equation (5) of Appendix III by the power ratios of column 8 of 
Table II. The values of Y, are then used in conjunction with the 
matrix of Appendix III to obtain the first approximation to the X, . 

In Fig. 24 the measured values of Y , expressed as decibel 
ratios with respect to beam 32, sre plotted as dots. A faired 


curve connecting the first approximation values of X, is indicated 


af 





Beam 


№ 


nN ኣጃ ድራ AM 


-J 


21 
22 


25 


TABLE I 


21 Ne 54 50.51 51 101.86 90150 51750 30 $0, 51,301.51, Gi ዝኽ 
43 49 49 48 k8 До э До 48 ho ДВ 48 49 49 48 46 47 ko Ag 
47 46 46 A6 47 46 Ат 48 47 45 46 47 AT 46 46 46 46 47 47 
49 45 48 48 49 48 48 ho ho 48 49 49 ho 49 ho 48 49 ho ho 
50 50 49 49 50 49 50 49 49 49 ho 49 50 50 49 49 50 51 49 
k7 k7 k7 47 k7 46 47 46 45 k7 k7 46 46 57 46 46 ki 46 ki 
49 50 49 49 49 48 49 50 49 54 51 50 50 49 50 49 50 50 50 
50 51 50 50 49 50 50 41 50 50 50 50 51 49 50 50 51 50 50 
48 49 49 48 lo 49 50 ho Ж 49 49 49 49 48 48 47 ho 50 48 
aO S 0750 50151 52.50 5) BOL5T BO poo SASA То 
AG 48 47 46 47 46 48 47 48 46 k7 k7 48 48 k7 k7 k8 48 46 
b2 44 43 44 45 44 45 45 45 43 44 44 45 44 44 43 ki 45 k3 
k2 hh 45 45 43 42 44 45 hh 45 44 45 43 43 hx hz 44 45 АҢ 
41 43 43 43 43 До 44 45 43 41 42 43 42 k2 42 42 43 44 h3 
46 49 10 49 49 48 ho 50 48 47 49 ho 49 49 48 47 48 Lo 47 
50 53 54 54 54 51 53 54 53 50 54 51 52 52 52 52 51 52 52 
Ат 50 ho 49 ho 49 49 ho 48 56 52 50 50 49 50 49 49 50 49 
5436 51.56 56 94 57"56 56 56 576 ВО Dp ws po po Юя 
48 51 51 50 50 50 51 50 50 49 49 50 51 50 49 48 50 52 49 
50 52 51 50 51 50 51 5% 52 51°50 50 49 50 51 51 50 30 51 
5%. 56 51 55 57756 5555095950158 0690000790500 
50 54 51 51 51 50 52 53 54 5252 25098 52 о а 
64 64 64 65 64 63 63 65 64 61 66 65 66 64 60 65 65 64 65 
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Beam 
24 
25 
26 
27 
28 
29 


32 


TABLE I (CONT,) 


62 62 62 61 61 61 62 61 60 59 61 62 62 60 62 59 60 
ШІ 550 1 500 50 Ga P150 MP gn 5251 51 5172 
ምክ FO 96 5655.96.59 57 56958, 'ጋፓማ ማዕ 37857 
60 29 29 60 2ይ 56 55 25 55 96 58 58 50 99 38 57 56 
60 59 60 61 60 59 59 59 61 58 60 61 60 58 58 57 58 
66 62 62 63 62 6h 65 6h 63 62 65 66 62 62 64 60 61 
69 74 67 65 64 64 66 65 61 64 66 69 67 68 63 64 70 
64 63 62 64 61 64 62 66 62 60 62 64 61 65 62 62 61 
77 78 76 74 74 76 75 77 73 74 75 76 72 74 75 16 77 


63 59 
21 51 
60 57 
60 57 
62 58 
65 63 
63 63 
61 59 
(game 


19 





Corr 


Bean Aver Aver Aver* Angle re 32 


1 22. 52 
2 h; ho 
3 16 18 
4 40 51 
5 48 52 
6 ls 48 
т м Юю 
8 an» 20 
9 48 47 

10 4B 51 

11 45 4 

12 13 45 

1 2 5 

1 44 47 

15 à; + 

16 yh М 

17 22 E 

18 k7 Аё 

19 45 4 

20 he 48 

21 61 «2 

22 SE "8 

23 56 52 

24 52 52 

25 552 =- 

26 61 595 

27 54 58 

28 54 54 

29 — 

20 59 64 

21 27" 59 

32 67 67 


TABLE II 

Vert Elev Ratio Ratio 

re ave 
Eo 225 603 999 
Se 1429 549 909 
h7 116" 495  .821 
16 101.59 .475  ,784 
E^ — 78.5? .611- 19080 
55 "649 sm 1.122 
60 589 .797 1.320 
67 0° 1.000 1.655 


Ratio 2. 
Re 32(db) Re av.(db) 


መ 


-2.61 


-2.06 


245 


-1.74 


0 


Ratio 


0.0 


0.42 


-1.06 


0.05 


0.46 


1,221 


2,19 


Y. 


/ 


„0522 


40297 


«0276 


0248 


„0519 


„0576 


„0417 


„0556 


*Vertical Average - The averaged output of all beams of equal eleva- 


tion angle. 
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FIGURE 24 TRUE FIELD VERTICAL DISTRIBUTION 
HERATIVE SOLUTION 





да 


by dashed lines. This curve is used to select values of X, in 10 
degrees steps to complete the numerical integration of equation (1) 
of paragraph 6.1. The integration of equation (1) was accomplished 
with incremental steps of 10 degrees in elevation and azimuth. 
Figure 25 represents the results of all canputations with the ex- 
ception of the integration over the variable, 2. The ordinate is 
G(ZQ) X (O)sinZ, the abscissa,  . The 8 curves represent the 8 
vertical beams and is included to illustrate the fact that the beam 
response contributes significantly in all directions and emphasizes 
the importance to this process of knowing the beam patterns exactly. 
The curves of Fig. 25 are then integrated by planimeter and the 
decibel ratios between a given beam and the beam direeted at the 
surface obtained. The results are plotted as X sin Fig. 24. The 
decibel difference between dots and Y 5 is measured and twice that 
value applied as a correction to the first assumption, dashed curve, 
of Fig. 24. The dotted curve results. The numerical integration 
process is repeated yielding the open circles. A final iteration 


yields the solid curve which is the curve of Fig. 18. 








FIGURE 25. G(8)X(8) SINO VS 8 
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